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A Nonconforming Finite Element Approximation for the von 

K arm an Equations 
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Abstract 

In this paper, a nonconforming finite element method has been proposed and analyzed for the 
von Karman equations that describe bending of thin elastic plates. Optimal order error estimates 
in broken energy and H 1 norms are derived under minimal regularity assumptions. Numerical 
results that justify the theoretical results are presented. 
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1 Introduction 


Let fl C ffi 2 be a polygonal domain with boundary dfl. Consider the von Karman equations for the 
deflection of very thin elastic plates that are modeled by a non-linear system of fourth-order partial 
differential equations with two unknown functions defined by: for given / £ L 2 (f2), seek the vertical 
displacement u and the Airy stress function v such that 


A 2 u 

A 2 v 


= [u, v] + / 
= -\[u,u] 


in 12 


with clamped boundary conditions 


u = 


du 

dv 


dv 

= «= — = 0 on dfl, 
dv 


( 1 . 1 ) 


( 1 . 2 ) 


where the biharmonic operator A 2 and the von Karman bracket [■, •] are defined by 


A p :— Pxxxx “b ^ pxxyy “b ^Pyyyy &nd [?7,y] I— 'HxxXyy ~b VyyXxx ^VxyXxy — Cof(D 7/) D y, 


coi(D 2 ri) denotes the co-factor matrix of D 2 n and v denotes the unit outward normal to the boundary 

on of n. 


Depending on the thickness to length ratio, several plate models have been studied in literature; 
the most important ones being linear models like Kirchhoff and Reissner-Mindlin plates for thin and 
moderately thick plates respectively; and non-linear von Karman plate model for very thin plates. 
Many practical applications deal with the Kirchhoff model for thin plates in which the transverse 
shear deformation is negligible. On the other hand, the Reissner-Mindlin plate model for moderately 
thick plates takes into consideration the shear deformation. The displacements of very thin plates are 
so large that a non-linear model is essential to consider the membrane action. The assumptions made 
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in the von Karman model are similar to those of Kirchhoff model except for the linearization of the 
strain tensor, which in fact, leads to the non-linearity in the model. 

For the theoretical study as regards the existence of solutions, regularity and bifurcation phenomena 
of von Karman equations, see n sm m na and the references therein. Due to the importance of 
the problem in application areas, several numerical approaches have also been attempted in the past. 
The major challenges are posed by the non-linearity and the higher order nature of the equations. 
The convergence analysis and error bounds for conforming finite element methods are analyzed in |12j . 
The papers [22] [25] and [24] investigate and analyze the Hellan-Hermann-Miyoshi mixed finite element 
method and a stress-hybrid method, respectively for the von Karman equations. In these papers, 
the authors simultaneously approximate the unknown functions and their derivatives. The papers 
[33 [13 HU deal with the approximation and error bounds for isolated solutions, thereby not discussing 
the difficulties arising from the non-uniqueness of the solution and the bifurcation phenomena. 


Over the last few decades, the finite element methodology has developed in various directions. 
For higher-order problems, nonconforming methods and discontinuous Galerkin methods are gaining 
popularity as they have a clear advantage over conforming finite elements with respect to simplicity 
in implementation. In this paper, an attempt has been made to study the von Karman equations 
using nonconforming Morley finite elements. The Morley finite element method has been proposed 
and analyzed for the biharmonic equation in m and for the Monge-Ampere equation in [23] . In [26] , 
a two level additive Schwarz method for a non-linear biharmonic equation using Morley elements is 
discussed under the assumption of smallness of data. The C° interior penalty method, a variant of the 
discontinuous Galerkin method has been used to analyze the Monge-Ampere equation in [91. 


The solutions u, v of clamped von Karman equations defined on a polygonal domain belong to 
Hq(0) n I? 2+Q (fl)[6], where a £ (|, 1] referred to as the index of elliptic regularity is determined by 
the interior angles of 12. Note that when fl is convex, a = 1. This paper discusses a nonconforming 


finite element discretization of (1.1)-(1.2) and develops a priori error estimates for the displacement 
and Airy stress functions in polygonal domains with possible corner singularities. To highlight the 
contributions of this work, we have 


• obtained an approximation of an isolated solution pair (u,v) of (1.1)-(1.2) using nonconforming 
Morley elements; 

• developed optimal order error estimates in broken energy and H 1 norms under realistic regularity 
assumptions; 


• performed numerical experiments that justify the theoretical results. 


The advantages of the method are that the nonconforming Morley elements which are based on piece- 
wise quadratic polynomials are simpler to use and have lesser number of degrees of freedom in com¬ 
parison with the conforming Argyris finite elements with 21 degrees of freedom in a triangle or the 
Bogner-Fox-Schmit finite elements with 16 degrees of freedom in a rectangle. Moreover, the method 
is easier to implement than mixed/hybrid finite element methods. 

The difficulties due to non-conformity of the space increases the technicalities in the proofs of error 
estimates. Moreover, one loses the symmetry property with respect to all the variables in the discrete 
formulation for nonconforming case. An important aid in the proofs is a companion conforming 
operator, also known in the literature as the enriching operator which maps the elements in the 
nonconforming finite element space to that of the conforming space. Also, as proved in m for the 
biharmonic problem, it is true that when Morley finite elements are used for the von Karman equations, 
the L 2 error estimates cannot be further improved. This is evident from the results of the numerical 
experiments presented in Section 

The paper is organized as follows. Section [l] is introductory and Section [2] introduces the weak 
formulation for the problem. This is followed by description of nonconforming finite element formula¬ 
tion in Section [3 Section [4] deals with the existence of the discrete solution and the error estimates 








2 WEAK FORMULATION 


3 


in broken energy and H 1 norms. The results of the numerical experiments are presented in Section [ 5 J 
Conclusions and perspectives are discussed in Section [G] The analysis of a more generalized form of 
fnj)-([L2j) is dealt with in Appendix A. 

Throughout the paper, standard notations on Lebesgue and Sobolev spaces and their norms are 
employed. We denote the standard L 2 scalar or vector inner product by (•, •) and the standard norm 
on for s > 0 by || • || s . The positive constants C appearing in the inequalities denote generic 

constants which may depend on the domain Q but not on the mesh-size. 


2 Weak formulation 

The weak formulation corresponding to fi~T|)-([L2| is: given / £ L 2 (fl), find u,v £ V := Hq(0) such 
that 


a(u, <pi) + b{u,v, <pi) + &(u,u, ipi) = l(ifi) V^iGH (2.1a) 

a[v, ip 2 ) — b(u, u, ip 2 ) = 0 \/ip 2 eV (2-lb) 


where V x, <P £ V, 

a(v,x) ■= [ D 2 rj : D 2 X dx, b(r],x,<p) '■=\ [ coi(D 2 r 1 )Dx ■ Dip dx and l(ip) := 

Jn 1 Jn 

Note that 6(-,-) is derived using the divergence-free rows property T5, .231- Since the Hessian matrix 
D 2 r] is symmetric, coi(D 2 rj) is symmetric. Consequently, 6(-,-,-) is symmetric with respect to the 
second and third variables, that is, b(rj,^,ip) = b(rj,tp ,£). Moreover, since [-,-] is symmetric, &(-,•,-) is 
symmetric with respect to all the variables in the weak formulation. 

An equivalent vector form of the weak formulation which will be also used in the analysis is defined 
as: for F = (/, 0) with / G L 2 (fl), seek T = (u, v) G V := V xV such that 


A(^,$) +5(4',^,$) = L($) V$£V (2.2) 

where V S = (£i,6),@ = ( 61 , 82 ) and $ = (^ 1 ,^ 2 ) G V, 

A(0,d>) := a(0i,v3i) +a(d 2 ,ip 2 ), (2.3) 

5(S,0,$) :=b(it 1 ,6 2 ,ipi)+b(£ y2 ,0i,Vi)-b(£,i,9i,ip 2 ) and (2.4) 

T($) := (/,¥>,). (2.5) 

It is easy to verify that the bilinear forms A(-, •) and !?(-, -, •) satisfy the following continuity and 
coercivity properties. That is, there exist constants C such that 

A(0,$) < C10||| 2 I $ 1 2 V0, $ G V, (2.6) 

A(0,0) > C 10||| 2 V0GV, (2.7) 

H(5,0,4>) < ClHla |||0| 2 !$|| 2 VS, 0, $ G V, (2.8) 


where the product norm ||| < I >|||2 := \J A(4>, <h) V$ G V. In the sequel, the product norm defined on 
(H s (O)) 2 and (L 2 (Ll )) 2 are denoted by |||-1|| s and |-||, respectively. 

For the results on existence of solution of the weak formulation, we refer to mmmm- M ° re 
precisely, the weak solution T = (u, v ) of can be characterized as the solution of the operator 

equation /'F = T’l' defined on V where T is a compact operator on V and / is an identity operator on 
V. In p], it has been proved that there exists at least one solution of the operator equation. Also, 
the uniqueness of solution under the assumption on smallness of the data function / has been derived. 
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Figure 1: Morley element 


In this paper, we follow T2j and assume that the solution 4/ = (it, v) is isolated. That is, the 
linearized problem defined by: for given G = ( 91 , 32 ) £ (T 2 (fl)) 2 C V', seek 0 = (9 1 , 62 ) £ V such that 

.4(0,$) = (G,$) V $ £ V (2.9) 

where .4(0, $) := 4(0, $) + U('k, 0, $) + B(0, T, <!>) is well posed and satisfies the a priori bounds 

|||©||| 2 <C|||G|||, |||©||| 2+a < C |||G||| (2.10) 

where a is the index of elliptic regularity. 

3 Nonconforming Finite Element Method (NCFEM) 

In the first subsection, the Morley element is defined and some preliminaries are introduced. In 
the second subsection, nonconforming finite element formulation for von Karman equations and the 
corresponding linearized problem are presented. Some properties and auxiliary results necessary for 
the analysis are discussed in the third subsection. 

3.1 The Morley Element 

Let Th be a regular, quasi-uniform triangulation [TQ S [13] of O into closed triangles. Set hr = 
diam(T) VT £ Th and h = max^g^ hr■ For T £ Th with vertices = ( Xi,yt ), i = 1,2,3, let 
TO 4 ,TO 5 and me denote the midpoints of the edges opposite to the vertices 01,02 and 03 respectively 
(see Figure [TJ. We denote the set of vertices (resp. edges) of Th by TI h (resp. fibj. For e £ <Eh, let 
h e = diam(e). 

Definition 3.1. [13] The Morley finite element is a triplet (T, Pt, $t) where 

• T is a triangle 

• P T = P 2 (T) is the space of all quadratic polynomials on T and 

• $ T = 1 are the degrees of freedom defined by: 

dv 

<t>i(v) = v(a,i), i = 1,2,3 and = —(m,), i = 4,5,6. 

The nonconforming Morley element space associated with the triangulation Th is defined by 

Vh := £ L 2 (VL) : p\t £ PiiT) VT £ Th, p is continuous at the vertices {ai}i =1 of the triangle 

and the normal derivatives of ip at the midpoint of the edges {mi }® =4 are continuous, 

ip = 0 at the vertices on dLl, ^ = 0 at the midpoint of the edges on dfl\. 

ov ) 
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For <p £ Vh and $ = (<pi, ip 2 ) £ Vh '■= Vh x Vh, the mesh dependent semi-norms which are equivalent 
to the norms denoted as \<p\ 2 ,h and |||$||| 2 h , respectively, are defined by: 

W\\,h := E Mi,Ti Wl 2 ,/i := IV^l 

TeT h 

Also, for a non-negative integer m, 1 < p < 00 and ip £ W m ' p (D\ Th) 

M m,p,h := E M m,p,T> IMIm.pA := E IMHm,p,T> 

TeTh TeTh 

and for p = 00 

|^M,oo,/i ■— max |^|m,oo,T5 IMU.oo.h •— max ||(p||m,oo,T) 

TeTh TeTh 

where | • | m ,p,T and || • || m ,p,T denote the usual semi-norm and norm in the Banach space W m ’ p {T) 
and W m,p (D;Th) denotes the broken Sobolev space with respect to the mesh Th- For <!> = (<pi,p 2 ) 
with <pi,<p 2 £ W m ’ p (n-T h ), define |||$|||^ >p>/l := \p>i\ 2 m ^ h + |<p 2 | 2 m , p , h - When p = 2, the notation is 
abbreviated as | • | mj h and || • || 

m,h • 


3.2 Nonconforming Finite Element Formulation 


The NCFEM formulation corresponding to (2.la)-(2.1b) can be stated as: for / £ L 2 (fl), seek 
( Uh,Vh ) € Vh such that 


a h {u h ,<pi) + b h (u h ,v h ,(pi) + b h (v h ,u h ,(pi ) = l h (<p 1 ) V<pi £ V h 
ah{v h -,<p 2 ) - b h (u h ,u h ,(p 2 ) = 0 

where V 77 , %, ip G V h , 


(3.1a) 

(3.1b) 


ah(v,x) '■= V / D2 V-D 2 xdx, b h (ri,X,v)'■= t E / cof (D 2 r])Dx- Dip dx and 
Ten Jt 2 — Jt 

h(<P) ■= E dx - 

Jt 


TeTh' 


TeTh' 


As in the continuous formulation, the discrete form bh(-,-,-) is symmetric with respect to the second 
and third variables. However, unlike in the conforming case HU, b h (; y) is not symmetric with 
respect to the first and second variables or the first and third variables. The equivalent vector form 
corresponding to (3.la)-(3.1b) is given by: seek = (uh,Vh) G Vh such that 


A h (^ h ,<S>) + Bhi’Fh^h,^) = L h (&) V$GV k 
where V E = (£ 1 , 6 ), 0 = (0i,0 2 ) and $ = ((pi,<p 2 ) £ Vh, 

A h (B,®) := 0^(01,<pi) + a h ( 9 2 l ip 2 ), 

B h (E, 0,<F) := bh{£i,9 2 ,<pi) + b h (t 2 ,0 1 ,ip 1 ) - b h {£i,9 lt (p 2 ) and 

L h(§) ■= E / fVi dx. 

Ten 


(3.2) 

(3.3) 

(3.4) 

(3.5) 


The nonconforming finite element formulation corresponding to (2.9) reads as: for given G G 
(L 2 (fl)) 2 , find Oh £ Vh such that 


A(0„,$) = (G,$) V $ £ V ;i 


(3.6) 


where A h (O h , $) := A h {O h , $) + Bh{F-. @ h , $) + B h (O h , 'F $) and A h (-, ■), B h (-, -, ■) are defined in 
(3.3) and ( |3.4| ), respectively. 
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3.3 Auxiliary Results 

In this subsection, some auxiliary results which are essential for the analysis are stated. 

Lemma 3.1. (Integral average) E The projection P e : L 2 (T) —> Po(e) defined by P e p = — j p ds, 

J e 

satisfies 

\W - P e <p\\ 0 ,e < Ch¥ 2 \<p\ 1>T Vp£H 1 (T). (3.7) 

Lemma 3.2. (Intervolant) ITU [W\, I T O) Let : V —> 14, be the Morley interpolation operator defined 
by: 


(lift A) (p) = p(p) 



VpS^Jft, 

dp 

ds Vee£ft. 


Then for p £ H 2+a (0), a € (0,1], it holds: 


||V nft^||m,p,/i 


< Ch 1+a ~ m+ p 


IMk+cn 


0 < m < 2,1 < p < oo. 


For simplicity of notation, the interpolant of <I> £ V is denoted by IIft$ and belongs to Vft. 

Lemma 3.3. (Enrichment function) 111) Let V c be chosen as Hsieh-Clough-Tocher macro element 
space num which is a conforming relative of the Morley finite element space Vft. For any p £ Vft, 
there exists EhP £ V c C V such that 

(h T 4 |</> — Ehp\o >T + h T 2 \p — Ehp\\^ T ) + l-Eft^l^ft < CM 2 ,ft- (3-8) 

T&T h 

Again, for $ £ Vft, the enrichment function corresponding to $ denoted by Pft$, belongs to V. 

In the next lemma, we establish an imbedding result. A similar result has been proved in [261 
Lemma 3.1] for the case of convex polygonal domains. However, for the sake of completeness, we 
provide a detailed proof for the case of polygonal domains. Note that only the edge estimation in 
is different from the proof in [26] . 

Lemma 3.4. (An imbedding result) For p £ Vft, it holds: 


(3.12 


Ml, 4, ft < CM2,ft. 


Proof. The tangential and normal derivative of p £ Vft are continuous at the midpoint of each edges 
of T £ 7ft. That is p x , p y £ Sh where Sh is the nonconforming Crouzeix-Raviart finite element space 
defined by 


Sh := £ L 2 (0) : w\t £ -Pi(T) VT £ 7ft, w is continuous at the midpoints of 

the triangle edges and w = 0 at the midpoint of the edges on 9fl|. 

It is enough to prove |ie|o, 4 ,ft < M|i,ft Vw £ Sh- 

Consider the auxiliary problem: given 9 £ seek £ such that 

= 0 in 12, ( = 0 on dLl. (3-9) 

The solution satisfies the following a priori bounds 


liei| 1 <c'||0||_ 1 , M| 1+7 <c|M, 


( 3 . 10 ) 
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where 7 £ (^,1] denotes the elliptic regularity of the problem (3.9). Let //,,£ £ <S\ be an interpolant 
which satisfies the estimate mm 


ic - ht\o,h + h \t, - uh,h < c'/ i i +^ieiii +7 . 


A multiplication of (3.9) with w and a use of Green’s formula leads to 

(0,w) = (-A£,w) = ^(V£,Vw)- j f ds 


TeTh 


TeTh' 


The boundary term can be estimated as follows: 

dt; 


E / ds = E E I ~ Pe^ ds. 


TeT h 


idT dv 


Ter h ecar- 


Si' 


r d 

Since I (w — P e w ) ds = 0 Ve £ Gp, and „ I hi, is a constant over each edge, we obtain 
Je OU 

E E ! ^{w-PeW) ds= Y E [ i-M - Ih 0(w - p e w) ds 


TeTh e<ZdT ' 


TeTh eCdT- 


< E E ll^- / ^l|l,e|k-- P e^l|o,e- 

TeT h eCOT 


A use of trace theorem, Lemma |3.1| and (|3.11|) leads to the estimate 

<ch^u\\i+Mhh- 


E 

T&Th' 




dT 


dv 


ds 


Therefore, the a priori bounds in (3.10) yields 


(0,w) < (|e|i + C^||e|| 1+7 )|«;| lih < C'(|| 6»||_ 1 + h?\\0\\)\w\ hh . 


A choice of 6 = w 3 in (3.14) leads to 

Ho, 4 l / l <G(lk 3 ll-i + ^ll^ 3 ll)kliA- 

A use of inverse inequality yields 

ll w3 || = IMI 3 6(n) < Ch-i\\w\\ 3 LHny 

Also, Holder’s inequality and the imbedding result L 4 (fi) lead to 

(w 3 ,0 = IMliWI^IU 4 ^) ^ l^li => ll™ 3 ||-i < C'lMli^n)- 

Hence, a use of (3.16) and ( 3.17[ ) in ( |3.15[ ) leads to the required result 

Mo,4 ,h < C( 1 + /l 7 _ 5) \w\i ih < C\w\ lt h- 

The next lemma follows from El Lemmas 4.2 & 4.3]. 

Lemma 3.5. (Bounds for Ah{-,-)) (i) Let x G (H 2+a (0)) 2 and <f> £ Vh■ Then, it holds 

A h {Xi E h <& — $) < Ch a ||x||| 2+a m*,h- 
(ii) Further, for x G (H 2+a (0)) 2 and $ £ (Hq(0)) 2 n (H 2+a (0)) 2 , it holds 


A h { X ,Hh^-^)<Ch 2 


III $||| 

2+a Ill^lll2+a * 


(3.11) 


(3.12) 


(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 


□ 
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A use of the definition of generalized Holder’s inequality and Lemma 3.4 leads to a bound 

given by 

Bh{ S, 0, $) < C b |S| 2ih |||0||| 2jh m 2 , h , (3-18) 

where C b is a positive constant independent of h. 

Lemma 3.6. (A bound for B b (-, •, •)) For 2 G (H 2+a (0)) 2 and @,$€ V + Vh, there holds 
Bh{ S,0,S) < CI|S|| 2+a ||0| li4 , fc ||$|| lih < C|||2||| 2+a |||0||| 2i , ||$| lih . 

Proof. Consider 

bh.(v,X,<P) = o / ((VyyXx - VxyXy)Tx + (VxxXy ~ VxyXx)p v ) dx. (3.19) 

2 TeT h Jt 

For r] G H 2 +a (D), a use of generalized Holder’s inequality and the imbedding result H 2 +a (Q) c —^ 
fU 2 ’ 4 (fl) leads to an estimate of the first term on the right hand side of (3.19) as 


hyyXxTx dx < ( \t]\ 2 ,4,t) ( l^b 


2 ,T 


TGTh ' 


\TdTh 


\TdT h 


\TdT h 


< |M|w 2 . 4 (n) |x|l,4,h Ml ,h < CIMh+a |x|l,4,/i \p\l,h- 


Similar bounds hold true for the remaining three terms in (3.19). Hence the required result follows 
using the definition of B b {-, ■, •) and Lemma 3.4 □ 


Remark 3.1. Using a proof similar to that of Lemma 3.6 it can be deduced that for 2 G (H 2 +a {0 )) 2 
and 0, $ G V + Vh, there holds 

B h { 2,0,$) <cm 2+a I© 11^ mx Ah < c , isi 2+a 101 m \ 2 , h . 


Using the definition of bh {•), an integration by parts and a use of (3.7), the following lemma 
holds true. 

Lemma 3.7. (An intermediate result) For rj G (V D H 2 +a (D)) + Vh and x, V £ Vh, it holds 
h{v, X, = b h {x, V, V) + o Y] / (VxXy - VyXx) Vip-T ds 

2 TeTfe •* dT 

where r is the unit tangent to the boundary dT of the triangle T. Moreover, 

V 77 , X) ^ £ V h , 7 : ^ f (VxXy - VyXx) Vif ■ t ds < Ch\ri\ 2 , h \x\i,oo,h M 2 ,ft- (3-20) 

Z J dT 


TeT h 


Remark 3.2. A use of Lemma [3.7[ Remark 3.1 and imbedding result H 2 +a (Il) ► W 1 , 0 O (fi) leads to: 
for Z h ,$ h G V h and £ G {H 2 +a (0)) 2 , 


I B h (E h ,£, < h/i)| < I B h (£, S h , $h)| + Ch 12 ft | 


2 ,h 1 ^ 12 +o; lll^h| 2 ,h • 


(3.21) 


The next lemma which will be used to establish the well posedness of the linearized problem (3.6), 
follows easily under the assumption that 4/ is an isolated solution of ( 2 . 2 ). 
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Lemma 3.8. (Well posedness of dual problem) If ’L is an isolated solution of (2.2), then the dual 
problem defined by: given Q £ (i(f _ 1 (J2)) 2 , find C £ V such that 


A(<RC) = (<2, < ^ , ) V$ev 

is well posed and satisfies the a priori bounds: 


1 2 <g 


l 2+ « < C 


(3.22) 

(3.23) 


(Q,<P) 


where a denotes the elliptic regularity index and |< 3 |||_i : = sup 

Q))= 

Since the Morley finite element space Vh is not a subspace of V and the discrete form b h {-,-,-) 
is non-symmetric with respect to first and second or first and third variables, we encounter addi¬ 
tional difficulties in establishing the well posedness of the discrete problem (3.6) in comparison to the 
conforming case. 


Theorem 3.9. (Well posedness of discrete linearized problem) If'F is an isolated solution of (2.2), 
then for sufficiently small h, the discrete linearized problem (|3.6|) is well-posed. 


Proof. The space Vh being finite dimensional, uniqueness of solution of (3.6) implies existence of 
solution. Uniqueness follows if an a priori bound for the solution of (3.6) can be established. That is, 
we aim to prove that 

< cr |||G]|| (3.24) 


for sufficiently small h. For <f> £ Vh, using Lemma 3.6 and Remark 3.2 the following Garding’s type 
inequality holds true: 

A h {®, $) = A h (<S>, $) + B h OF, <R <f>) + B h (<S>, R $) 

> mlh -G|*||| 2+a ll^l 2 ,h 


Substitute <f> = 0^ in (3.6) and use (3.25) to obtain 


Note that 


m hh -ch i$i 2ih i'Fi 2+Q m 2 , h - 

(3.25) 

: + Ill’ll 2+0 111,*. + G1). 

(3.26) 

'ft- IIIi,/i — Gh 0ft 2 ,ft + P?ft0ft k • 

(3.27) 


Now we estimate |||IR0/,||| 1 . Choose Q = —A E h @h and 4) = E h @h hr (3.22) and use (3.6 1 to obtain 

l\E h <d h i\l = A(E h O h , C) = A h (E h O h , C - IRC) + A h (E h Q h , IRC) 

= A h (E h e h , C - IRC) + Ah(E h O h - Oh, IRC) + (G, IRC) 

= A h (E h o h — o h , C) + A h (o h , C - IRC) + B h (F, E h 0 h , C - IRC) + B h (E h o h , R C - IRC) 
+ B h (V, E h O h - Oh, IRC) + B h {E h Oh - 0 h , R IRC) + (G, IRC) (3.28) 

A use of Lemmas 3.2 3.3 3.5 ( 3.18| ), (3.23), Remarks 3.1 and |3.2| leads to 

1^0,11? < C (V ffl0 h || 2ih ICIIa+a + h ° I^III 2 I0J.II2,/, HClIa+a + h |||R|| 2+Q ||0 fc | 2 , h |||CIII 2 + |||G|| |||C||| 2 


< c [h 


Therefore, 


(h a |||0/J 2) / i + III G III) < C ( h a HI ©ft. 

\\E h e h tu h <c(h a ||0 h |o h + 


12 ,h 


i\E h O h |r. 


Now, (3.26)-(3.29) yield 

l|0h|| 2 ,/, < C.h a |||©h||| 2) / l + C IGI. 

That is, 10/j .||| 2 h < C |||G||| for a choice of h < h\ = ( 5 ^)° with a £ (|, 1]. 


(3.29) 


□ 
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Remark 3.3. If is an isolated solution of (2.21, then for sufficiently small h, the discrete linearized 
dual problem: given Q G (L 2 (H)) 2 , find G Vh such that 


= ($,$) vsgv* 

is well posed. The proof is similar to that of Theorem |3.9|and hence is skipped. 


(3.30) 


4 Existence, Uniqueness and Error Estimates 

In view of Theorem |3.9| and Remark 1 3.3 1 the bilinear form Ah{-, •) : Vh x Vh —■> R defined by 

-U(0, $) = A h { 0, $) + , 0, $) + B h (Q, U $) (4.1) 

is nonsingular on Vh x Vh- 

The next lemma establishes that the perturbed bilinear form A;, (•,•), constructed using IR’I' is 
also nonsingular. Though a similar result is proved in m for the conforming case, we provide a proof 
here for the sake of completeness. 

Lemma 4.1. (Nonsingularity of perturbed bilinear form) LetUh 4' be the interpolation of'It as defined 
in Lemma \3.2\ Then, for sufficiently small h, the perturbed bilinear form defined by 

Ah{ 0, <b) = Ah(0, $) + , 0, $) + B h (Q, n^, $) (4.2) 

is nonsingular on Vh xVj,, if ( |4.1[ ) is nonsingular on Vh x Vh- 

Proof The bilinear form Ah ■ Vh x Vh —> R is bounded and satisfies 

SUp A(©,$) > /3|M 2 h> SUp A(0,$)>/8|0|||2/ lI 

l|e| a ,h=l ’ l*la,ii=l 


(see H3.18D) yields 


where /? > 0 is a constant. For 'PgV + Vj, a use of the above properties of Ah{-, ■) and continuity of 
fields 

sup A h (Q, $) + B h ()I> - T, 0, $) + B h (Q, V - $) 

M a . h =i 

sup A(0,4 > )- sup (s ft (§, 0, $) + .B fc (0, §,$)') 

fill , >,=1 111 * 111 , ,=1 V 7 


> 


provided 




— HI 0 ft. HI 2 ,*. 
P 


4- 


2,/i 


2 ,fc > 2 Il e ll 2 ,fc > 


— 77— ■ Such a choice is justified for sufficiently small h < hi (say), by setting 
2 ,h 4G 5 


’F = ’F — II/j7k and using Lemma 
required result. 


3.2 


Similarly, sup .4/ l (0,4>) > — |$1 2 h V <F G V/t. Hence the 


|e|| a .h=i 


□ 


4.1 Existence and Local Uniqueness Results 

Consider the nonlinear operator p : Vh —► Vh defined by 

Ah(ffiO), $) = L h ($) + , 0, $) + B h (Q, n h V, $) - B h {Q, 0, $) V $ G V h . (4.3) 

A use of Lemma |4.1| leads to the fact that the mapping p is well-defined and continuous. Also, any 
fixed point of p is a solution of (3.2) and vice-versa. Hence, in order to show the existence of a 
solution to ( |3.2| ), we will prove that the mapping p has a fixed point. As a first step to this, define 
Bfl(n h tt) := {$ G V h : ||$ - n h *1 2>h < i?}. 
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Theorem 4.2. (Mapping of ball to ball) For a sufficiently small choice of h, there exists a positive 
constant R(h) such that for any 0 € Vh, 

lie - M| 2i „ < R(h) =► !m( 0 ) - Mi 2ih < R W- 

That is, /i maps the ball B^(; 1 /)(n/ 1 'F) to itself. 


Proof. Since the bilinear form Ah {■, ■) is nonsingular, from Lemma 4.1 there exists $ € Vh such that 


p> L , = 1 and 

II 1112 , h 


§ IH©) - n„*||| 2ih < Ah(^Q) - M, $) 


Let Eh$ be an enrichment of $ (see Lemma 3.3). A use of (4.2), (4.3) and (2.2) yields 

A h (p(o) - m, $) = A(m(0),*) - M n ^, *) 

= M) + B h {U h A, 0, $) + B h (Q, M, $) - B h (0,0, $) - a„(m, $) - 2B h (M, m, $) 

= M - E,A) + (A h {*, M) - A h ( m, $)) + (B h ($, 4-, M) - B h (n h y, m, $)) 

+ I^(M - ©, 0 - M, #) =: Ti + T 2 + T 3 + r 4 . (4.4) 


Now we estimate {Tj}| =1 . Ti can be estimated using Lemma 3.3 and the continuity of Lh- Using 
Lemma |3.5| continuity of Ah(-, •) and Lemma |3.2[ we obtain 


T 2 < \A h (*,E h $)- A h (IL h *,$)\ < \A h ('V,E h $-$)\ + \A h ('V-Il h 'V,$)\<Ch a |$| 


2+a ' 


A use of Lemmas 3.6 3.3113.21 and (3.18) leads to 


T 3 < |5 /t (4s4sM)-^(M,M,4>)| 

< I B h {*, 'L,M - $) - B h (M - 4-,M, $) - Bfc(w ,m -*, #)| 


< C , fc|||$||| 2+Q I^1 2 |||M| 

<Ch a |tt|| 2 I|*||| 2+a . 


u+ch-mUa ii^iia winkh 


Finally, X 4 is estimated using (3.18) as 


T a < |B h (M - 0, 0 - U h A, $)| < C |||0 - Ml 


2 ,h ' 


A substitution of the estimates derived for Ti,T 2 ,T 3 and T 4 in (4.4) and an appropriate grouping of 
the terms yields 


IIH©) - < C x (h° 


'-Ml 


2 ,h 


(4.5) 


for some positive constants C\ independent of h but dependent on ||| v f'||| 2 +Q ,- A choice of h < / 13 , where 
h 3 = ( 4^2 ) ° i yields 4 Cfh a < 1 . Since |||0 — n /, 4'||| 2 h < R(h), for h < h 3 , a choice of R(h) := 2C\h a 


leads to 


M&) - Mill, <h < C x h a (1 + 4 (%h a ) < R(h) 


This completes the proof. 


□ 


Theorem 4.3. (Existence) For sufficiently small h, there exists a solution 'Vh of the discrete problem 
(3.2) that satisfies |4'/ ( — II/ l 4 , || 2 h < R(h), for some positive constant R(h) depending on h. 
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Proof. Lemma 4.2 leads to the fact that /i maps the ball (11^4') to itself. Therefore, an application 
of Schauder fixed point theorem [lH] yields that the mapping /i has a hxed point, say T i,. Hence, 4/^ 
is an approximate solution of (|3.2|) which satisfies ||| 'S/h ~ II/ 1 ' , QE r ||| 2 h < R(h). □ 


Theorem 4.4. (Contraction result) For ©i ,©2 G (11^41) with R(h) as defined in Theorem f.2 
the following contraction result holds true: 


IM0i) - M(02)||| 2ife < Ch a 10i - 0 2 || 2ifc , (4.6) 

for some positive constant C independent of h. 

Proof. For 0 1 ,0 2 G let /i(©i),i = 1,2 be the solutions of: 


A(/i(©,), $) = L h ($) + B h ( n h $, ©i, <S>) + B h (Q it n^, 4>) - B h (Q z , ©i, $) V 4> G V h . 


(4.7) 


The nonsingularity of Ah(-, •) yields a $ with 4> = 1. With (4.7) and (3.18), we obtain 


{ IIH©i) - m(© 2 )ll 2 ,/i < A(M©i) - m(© 2 ), *) 

= B h {U h V, ©1 - ©2, $) + B h (Q 1 - ©2, n h *, $) + 5^(02, ©2, $) - B fc (©i, ©!, #) 
= B h (0 2 - ©i, ©1 - n ;i T, $) + Bh (0 2 - n h *, ©2 - ©i, $) 

< c 1 ©2 - ©i| 2ih (|||01 - n ^|| 2 h + 1||02 - n h n 2 , h ) ■ 


Since ©i, ©2 G B^^n^Hi), for a choice of R(h) as in the proof of Theorem 
h, we obtain 

!M©i) - M©2)lll 2,h<Ch a |||02-0l| 2lh) 

for some positive constant C independent of h. This completes the proof. 


4.2 


for sufficiently small 

(4.8) 

□ 


Remark 4.1. (Local uniqueness) Let 4/ be an isolated solution of (2.2). For sufficiently small choice of 
h , Theorem 4.4 establishes the local uniqueness of the solution of (3.2). 


4.2 Error Estimates 

In this subsection, the error estimates in the broken energy and H 1 norms are established. 


Theorem 4.5. (Energy norm estimate) Let 4> and ’F/i be the solutions of (2.2) and (3.21 respectively. 

(4.9) 


Under the assumption that 4/ is an isolated solution, for sufficiently small h, it holds 


l*-*hl 2th <Ch° 


where a G ( 2 , 1] is the index of elliptic regularity. 
Proof. A use of triangle inequality yields 


ll*-*ft!2,h < Ill^-n^llU 




For sufficiently small h, Theorem |4.3| leads to 


HIM - < Ch a . 


Now, Lemma 3.2 , (4.11) and (4.10) establish the required estimate. 


(4.10) 

(4.11) 
□ 

















4 EXISTENCE, UNIQUENESS AND ERROR ESTIMATES 


13 


Theorem 4.6. (H 1 estimate) Let 4' and 44 be the solutions of (2.2) and (3.2) respectively. Assume 
that 4/ is an isolated solution. Then, for sufficiently small h, it holds 


!*-*/*! \i, h <Ch 2a , 


(4.12) 


where a £ (|, 1 ] is the index of elliptic regularity. 
Proof. A use of triangle inequality yields 


14/ — 4/ 


h\\\lh — 


< ii$ - u h m 


l,h 


HIM - $ 


Mill,ft ^ 


< ii$-n fc $i 


1 ,h 


Wp-Ehphb + WEbpi,, (4.13) 


where p = 11^4' — 44- A choice of Q = —AE^p and 4> = EhP in the dual problem (3.22) and a use of 
(| 2 . 2 |), ( 3.2h leads to 


(yE h p,\7E h p) = A h (E h pX) = A h {E h p- pX) + A h {p, <) 

= Ah{EhP — p, C) + Bh(X j EhP — p, C) + Bh(EhP — P , 4q £) 

+ A h (n„$ - $,C) + A h ($ - 44, c - 140 + A h (%n h c - 0 + MC - n h o 
+ (B h ( $, 14$ - $ h) C) + s fc (n h $ - $, o - B h ($, $, C) + B h ($ h) 44, 14c)) 

8 

=:E T - ( 4 - 14 ) 

i=l 


Xi is estimated using Lemma |3.5| and (4.11). T 4 and Tq are estimated using Lemma 3.5 T 5 is 
estimated using continuity of Aft(-,-), Lemma 3.2 and Theorem 4.5 The term T 7 is estimated using 
continuity of Lft and Lemma 3.2 T 2 is estimated using Remark 3.1 Lemma 3.3 and (4.11) as 


X 2 < 11^(4/, E h p - p, C)| < C ||$|| 2+a II E h p - pi, 


< Ch 1+a I $ I 


2+a I 


(4.15) 


T 3 is estimated using Remark 3.2 Lemma |3.3| (4.15) and (4.11) as 
T 3 < \B h (E h p-p,*X)\ < \B h (*,E h p-pX)\+Ch \\E h p - pl 2 h ||$| a+a | 


| 2 < Ch 1+ ° 


II $ I 


2 + 0 ! 


112 • 


Finally, a use of Remar ks |3.1| , |3.2[ Lemmas |3.2[|3.6[ Theorem 4.5 and (3.18) yields an estimate for 7g 
as 


t 8 = S h ($, 14 4> - 44, 0 + B h (n h ^ - $, 0 - s h { 4 -, 4>, c) + B h ( 44 , 44, 14 c) 

= B h (y, 144 - - $, 0 + s h (n h $ - 4-, 4>, c) 

+ B h ($, $ - 44, 0 + B h {i> - 44, $, C) - B h ($, $, 0 + 44,I 4 <) 

= B h (y, 14 4 > - 4 -, 0 + B h ( 14$ - 4-, 4 -, c) 

+ B h ($ - 44, $ - 44, 0 + B h { 44 - 4 /, 44, i 4 < - 0 + B/,(®, $/., I 4 C - 0 
<C'^“(IIICIIl2 + IIIClll 2+ J- 


A combination of the estimates T) to T§ and a priori bounds (3.23) for the linearized dual problem 
yields 


C S7E h p,\7E h p ) < C7fc 2 “ I-A^plt, < Ch Za lEnp^ =► \iE hP \\ x < Ch Za . 


(4.16) 


A use of Lemmas 3.2 3.3| (4.11) and the last statement of (4.16) in (4.13) completes the proof. □ 


4.3 Convergence of the Newton’s Method 

In this subsection, we define a working procedure to find an approximation for the discrete solution 
4/^. The discrete solution 4 4 of (3.2) is characterized by the fixed point of (4.3). This depends on 
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the unknown II/,Ik and hence the approximate solution for (3.2 1 is computed using Newton’s method 
in implementation. The iterates of the Newton’s method are defined by 


A h (n, $)+o(*r, + B h(*h $ ) = ^(*r ,> $ )+ L ^) v$ev k . (4.17) 


Now we establish that these iterates in fact converge quadratically to the solution of (3.2). 


Theorem 4.7. (Convergence of Newton’s method) Let 4/ be an isolated solution of (2.2) and let 40, 
solve (3.2). There exists p > 0, independent of h, such that for any initial guess 4-2 which satisfies 

2h < P, ffl*S - *fc| 2 K < -3- holds true. That is, the iterates of the Newton’s method 


11*2 -ib 


defined in (4.171 are well defined and converge quadratically to 4'/,. 

there exists 6 > 0 such that for each Z h £ V h satisfying ||| Z h — II/,4 '||| 2 h 


Proof. From Lemma 
the form 


4.1 


<<5, 


A h (Q, 4>) + B h (Z h , 0, $) + 0(0, Z h , $) 


(4.18) 

is non singular in Vh x V/,. From (4.11), for sufficiently small h, 104- — 4 '/ l ||| 9 h < Ch a . Thus h can 
be chosen sufficiently small so that §04- — 4 -/,||| 2 h < Define 


p := mm 


5 /3 


2 ’ 16C b 


(4.19) 


where j3 and O are respectively the coercivity constant of _4/,(-, ■) and boundedness constant of 0(‘, v) 
(see (|3.6|)). Assume that the initial guess 4'° satisfies |||4>^ — 4-2 ||| 2 h < p. Then 


|n h * - *h|| 2 h < 1^4- - *hia,/, + ph - * 


2 || 2h <* 


Since (4.18) is nonsingular, the first iterate 4-2 of the Newton’s method in (4.17) is well defined for the 
initial guess 4-2- Using the nonsingularity of (4.18), there exists 4> £ Vh such that |||$||| = 1 which 


satishes 




2 .h 


< A h ( $) + 0(4-2, - 4o $) + 0(*i - *2, *)■ 


A use of (4.17), (3.2), (3.18) yields 


Hence, 


\n - ** 


Aki^i - 4o $) + 0(4-2, 4-2 - 4- h , *) + B h (*£ - 4//,, 4-°, 4>) 

= B h { *2, ., *) + L h {$) - A h (* h , $) - 0(4-2, * h , $) - 0(4-/,, 4-2, T) 

= 0(4-2,4-2, $) + 0(4-/,, */,, $) - 0(4-2,40 $) - 0(*fc, 4-2,4>) 

= 0(4-2 - 4-*., 4-2 - 40$) < O 14-2 - 4-/, 1112 

2 


\2,h ' 


I ^ 8Cb 
1 2,h ~ /3 


4-2-4-J 


2 h . Since |||*2 — 4-/,||| < p < we obtain 


(4.20) 




2, h 


< ^ III 4-2-4-/, 


I < P ~ 

\2,h - 2 ■ 


(4.21) 


Since |||4-^ — 4-/,||| < p, the form (4.18) is nonsingular for O = 4-2- Continuing the process, we 


obtain 


||4-i - 4- 


h\\\2,h ^ 


< — . 


Moreover, proceeding as in the proof of the estimate ( |4.20 1 , it can be shown that 

l4-" +1 - 4-J" ^ /m 111 ' T ' n - Tl 1112 


' h - ^h||| 2ih < (8Cf,//3) III *2 — 4-ft, III 2 ^ . 


(4.22) 


(4.23) 


This establishes that the Newton’s method converges quadratically to 40 This completes the proof. 

□ 
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Remark 4.2. (Local uniqueness) The local uniqueness of solution of (3.2) also follows from Theorem 4.7 


We observe that the definition of p in (4.191 does not depend on h. From Theorem 4.7 it is clear that 
for any initial guess 4'° which lies in the ball of radius p with center at 4^, the sequence generated by 
(4.171 will converge uniquely to 4v In particular, if we choose the initial guess = LEU/, then the 
sequence generated by the iterates of the Newton’s method will also converge to T/ t which shows the 
local uniqueness of the solution 4V 


5 Numerical Experiments 


In this section, two numerical experiments that justify the theoretical results are presented. The 
implementations have been carried out in MATLAB. The results illustrate the order of convergence 


obtained for the numerical solution of (1.1)-(1.2) computed using the Morley finite element scheme. 


For a detailed description of construction of basis functions for the Morley element, see Ming & Xu 


We implement the Newton’s method defined in (4.17) to solve the discrete problem (3.2). 


5.1 Example 1 

In the first example, we choose the right hand side load functions such that the exact solution is given 

by 


u(x, y ) = x 2 (l — x) 2 y 2 { 1 — y) 2 \ v(x , y) = sin 2 (7ra;) sin 2 (7r y) 

on the unit square. The initial triangulation is chosen as shown in Figure |2ja). In the uniform 
red-refinement process, each triangle T is divided into four similar triangles (I] as in Figure §b). 

Let the mesh parameter at the TV-th level be denoted by hfj and the computational error by ejv. 
The experimental order of convergence at the TV-th level is defined by 

a N := log(e N -i/e N )/log(h N _i/h N ) = log{e N - 1 /e N )/log{2). 




(a) (b) 

Figure 2: Initial triangulation 7o and its red refinement 71 


Tables HI and [2] show the errors and experimental convergence rates for the variables and ■ In 
Figures 3]|4 the convergence history of the errors in broken energy, H 1 and L 2 norms are illustrated. 
The computational order of convergences in broken H 2 , H 1 norms are quasi-optimal and verify the 
theoretical results obtained in Theorems 4.5 and |4.6| for a = 1. The order of convergence with respect 
to L 2 norm is sub-optimal justifying the results in m that using a lower order finite element method, 
the order of convergence in L 2 norm cannot be improved than that of the H 1 norm. 
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# unknowns 

1 u - U h \ 2 ,h 

Order 

\u - u h \ 1<h 

Order 

||« - U h \\ L 2 

Order 

25 

0.874685E-1 

- 

0.102155E-1 

- 

0.386068E-2 

- 

113 

0.405787E-1 

1.1080 

0.257318E-2 

1.9891 

0.919743E-3 

2.0695 

481 

0.209921E-1 

0.9508 

0.732470E-3 

1.8127 

0.248134E-3 

1.8901 

1985 

0.106209E-1 

0.9829 

0.191118E-3 

1.9383 

0.636227E-4 

1.9635 

8065 

0.532754E-2 

0.9953 

0.483404E-4 

1.9831 

0.160158E-4 

1.9900 

32513 

0.266595E-2 

0.9988 

0.121213E-4 

1.9956 

0.401107E-5 

1.9974 


Table 1: Errors and convergence rates of Uh in broken H 2 ,H 1 and L 2 norms 


^ unknowns 

\v - V h \ 2 ,h 

Order 

\v - V h \l,h 

Order 

||V - V h \\ L 2 

Order 

25 

19.245671 

- 

2.140613E-0 

- 

0.770876E-0 

- 

113 

9.5043699 

1.0178 

0.569979E-0 

1.9090 

0.177898E-0 

2.1154 

481 

5.0549209 

0.9109 

0.161737E-0 

1.8172 

0.482777E-1 

1.8816 

1985 

2.5758939 

0.9726 

0.421546E-1 

1.9398 

0.123930E-1 

1.9618 

8065 

1.2944929 

0.9926 

0.106618E-1 

1.9832 

0.312076E-2 

1.9895 

32513 

0.6480848 

0.9981 

0.267351E-2 

1.9956 

0.781643E-3 

1.9973 


Table 2: Errors and convergence rates of Vh in broken H 2 , H 1 and L 2 norms 


5.2 Example 2 

Consider the L-shaped domain Q = (—1, l ) 2 \ ([0,1) x (— 1 , 0]) (see Figure [5j. Choose the right hand 
functions such that the exact singular solution j! 6 | in polar coordinates is given by 

u(r, 0) = ( r 2 cos 2 0 — 1 ) 2 (r 2 sm 2 0 — 1 ) 2 r 1 +a g a ^(0)\ v(r , 6 ) = w(r, 9), 

where w := and a := 0.5444837367 is a non-characteristic root of sin 2 (aw) = a 2 sin 2 (o;) with 

9 a,ui( 0 ) = ^ 1 ^ sin ((a - l)w)--^-y sin ((a + l)w)^ ^ cos ((a - 1 ) 0 ) - cos ((a + 1 ) 0 )) 

—j- sin ((a — 1 ) 0 )-j- sin ((a + 1 ) 0 )^ ^cos ((a — l)w) — cos ((a + l)w)) . 

Tables [3] and [4] show the errors and experimental convergence rates for the variables Uh and Vh■ The 
domain being non-convex, we do not obtain linear and quadratic order of convergences in broken 
energy and H 1 norms for displacement and Airy stress functions. 



# unknowns 

I u - u h \ 2 ' h 

Order 

I u - U h \l,h 

Order 

|it ^ U h \\ L 2 

Order 

17 

29.209171 

- 

6.363539E-0 

- 

2.769499E-0 

- 

81 

14.130192 

1.0476 

1.682747E-0 

1.9190 

0.693436E-0 

1.9977 

353 

7.5651300 

0.9013 

0.491659E-0 

1.7750 

0.200814E-0 

1.7879 

1473 

3.9620126 

0.9331 

0.146551E-0 

1.7462 

0.583024E-1 

1.7842 

6017 

2.0841141 

0.9267 

0.487106E-1 

1.5891 

0.179703E-1 

1.6979 

24321 

1.1252534 

0.8891 

0.187772E-1 

1.3752 

0.613474E-2 

1.5505 


Table 3: Errors and the experimental convergence rates for Uh in broken H 2 ,H l and L 2 norms for 
L-shaped domain 
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Figure 3: Convergence history of displacement for Example 1 


# unknowns 

|w - V h \ 2 ,h 

Order 

\v - D(, l ,h 

Order 

\\V - V h \\ L 2 

Order 

17 

24.759835 

- 

4.932699E-0 

- 

2.069151E-0 

- 

81 

15.293270 

0.6951 

1.779132E-0 

1.4712 

0.727981E-0 

1.5070 

353 

7.8509322 

0.9619 

0.483823E-0 

1.8786 

0.199644E-0 

1.8664 

1473 

4.0531269 

0.9538 

0.137278E-0 

1.8173 

0.557622E-1 

1.8400 

6017 

2.1219988 

0.9336 

0.439086E-1 

1.6445 

0.165699E-1 

1.7507 

24321 

1.1421938 

0.8936 

0.165883E-1 

1.4043 

0.545066E-2 

1.6040 


Table 4: Errors and the experimental convergence rates for Vh in broken H 2 ,H l and L 2 norms for 
L-shaped domain 


6 Conclusions & Perspectives 

In this work, an attempt has been made to obtain approximate solutions for the clamped von Karrnan 
equations defined on polygonal domains using nonconforming Morley elements. Error estimates in 
broken energy and H 1 norms are established for sufficiently small discretization parameters. Numerical 
results that substantiate the theoretical results are obtained. A future area of interest would be 
derivation of reliable a posteriori error estimates that drive the adaptive mesh refinements. 
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Figure 4: Convergence history of Airy stress for Example 1 



Figure 5: L-shaped domain and its initial triangulation 


References 

[1] J. Alberty, C. Carstensen, and S. A. Funken, Remarks around 50 lines of Matlab: short finite 
element implementation, Numer. Algorithms 20 (1999), no. 2-3, 117 -137. 

[2] M. S. Berger, On von Karman equations and the buckling of a thin elastic plate, I the clamped 
plate , Comm. Pure Appl. Math. 20 (1967), 687-719. 

[3] _, Nonlinearity and functional analysis, Academic Press, 1977. 

[4] M. S. Berger and P. C. Fife, On von Karman equations and the buckling of a thin elastic plate, 
Bull. Amer. Math. Soc. 72 (1966), no. 6, 1006-1011. 

[5] _, Von Karman equations and the buckling of a thin elastic plate. II plate with general edge 

conditions, Connn. Pure Appl. Math. 21 (1968), 227-241. 

[6] H. Blum and R. Rannacher, On the boundary value problem of the biharmonic operator on domains 
with angular corners, Math. Methods Appl. Sci. 2 (1980), no. 4, 556581. 




REFERENCES 


19 


[7] D. Braess, Finite elements, theory, fast solvers, and applications in elasticity theory, 3rd ed., 
Cambridge, 2007. 

[8] S. C. Brenner, Forty years of the Crouzeix-Raviart element, Numer. Meth. for PDEs 31 (2015), 
367 -396. 

[9] S. C. Brenner, T. Gudi, M. Neilan, and L.-Y. Sung, C° penalty methods for the fully nonlinear 
Monge-Ampere equation, Math. Comp. 80 (2011), no. 276, 1979-1995. 

[10] S. C. Brenner and L. R. Scott, The mathematical theory of finite element methods, 3rd ed., 
Springer, 2007. 

[11] S. C. Brenner, L.-Y. Sung, H. Zhang, and Y. Zhang, A Morley finite element method for the 
displacement obstacle problem of clamped Kirchhoff plates, J. Comput. Appl. Math. 254 (2013), 
3142. 

[12] F. Brezzi, Finite element approximations of the von Karman equations, RAIRO Anal. Numer. 12 
(1978), no. 4, 303-312. 

[13] P. G. Ciarlet, The finite element method for elliptic problems, North-Holland, Amsterdam, 1978. 

[14] _, Mathematical elasticity: Theory of plates, vol. II, North-Holland, Amsterdam, 1997. 

[15] L. C. Evans, Partial differential equations, vol. 19, American Mathematical Society, 1998. 

[16] P. Grisvard, Singularities in boundary value problems, vol. RMA 22, Masson & Springer-Verlag, 
1992. 

[17] J. Hu and Z. C. Shi, The best L 2 norm error estimate of lower order finite element methods for 
the fourth order problem, J. Comput. Math. 30 (2012), no. 5, 449-460. 

[18] S. Kesavan, Topics in functional analysis and applications, New Age International Publishers, 
2008. 

[19] G. H. Knightly, An existence theorem for the von Karman equations, Arch. Ration. Mech. Anal. 
27 (1967), no. 3, 233-242. 

[20] P. Lascaux and P. Lesaint, Some nonconforming finite elements for the plate bending problem, 
RAIRO Anal. Numer. 9 (1975), no. Rl, 9-53. 

[21] W. Ming and J. Xu, The Morley element for fourth order elliptic equations in any dimensions, 
Numer. Math. 103 (2006), 155-169. 

[22] T. Miyoshi, A mixed finite element method for the solution of the von Karman equations, Numer. 
Math. 26 (1976), no. 3, 255-269. 

[23] M. Neilan, A nonconforming Morley finite element method for the fully nonlinear Monge-Ampere 
equaton, Numer. Math. 115 (2010), 371 394. 

[24] A. Quarteroni, Hybrid finite element methods for the von Karman equations, Calcolo 16 (1979), 
no. 3, 271-288. 

[25] L. Reinhart, On the numerical analysis of the von Karman equations: mixed finite element ap¬ 
proximation and continuation techniques, Numer. Math. 39 (1982), no. 3, 371-404. 

[26] X. Xu, S. H. Lui, and T. Rahaman, A two level additive Schwarz method for the Morley non- 
conforming element approximation of a nonlinear biharmonic equation, IMA J. Numer. Anal. 24 
(2004), 97-122. 



7 APPENDIX 


20 


7 Appendix 


We consider one of the variants of von Karman equations which is important in practical applications 
and give a brief sketch of the extension of the analysis. Consider the following form of von Karman 
equations: 

A 2 u = [«, v] — jj Au + / 


A 2 u =— \[u,u\ 


in Q 


with clamped boundary conditions 


du dv n an 

u = — = v = — = 0 on oil, 

ov ov 


(7.1) 


(7.2) 


where p is a real parameter known as the bifurcation parameter and D denotes the flexural rigidity of 
the plate. The weak formulation of (LU - \7.2\ reads as: given F = (/, 0), find df £ V such that 

4('F,$) + 5('F,4',$) + £('F,4>) = L(4>) v$ev (7.3) 


(7.4) 


where ^4(-,-), L(-) are defined in ( |2.3| )-( 2.5) respectively, and £(•,•) is defined as 

£(©,$) = [ X79i -X7ipi dx V0 = (9i,0 2 ) and $ = {tpu<f> 2 ) G V. 

D Jn 

The corresponding nonconforming finite element formulation is given by: find £ Vh such that 

A h (* h ,$)+B h (* h ,* h ,$) + <i: h (* h ,$) = L h {$) (7.5) 


where A h {-,-), L h (-) are defined in (3.3)-(3.5) respectively, and £h(-, •) is defined as 


£*(©,$) = ^2 J V6b ■ V(/2i dx V0 = (6»i,6» 2 ) and <f> = (^1,^2) £ Vh- (7.6) 


TdTh 


For the newly introduced bilinear form £(•, •), the following boundedness properties hold true: 

£(©,$)< (71||©||| 1 |||$||| 1 V0,$£V (7.7) 

£h(0h.*h)<C , |0 h | lih |||$/i| 1>h V© h ,$ h eV h . (7.8) 


For the modified problem (7.3), the linearized problem (see (3.61) is defined by: for given G £ (L 2 (fl)) 2 , 
find 0 £ V such that 

.4(0, 4>) = (G, d>) V$£V (7.9) 

where 

71(0, $) := 4(0, $) + 5(>F, 0, $) + 5(0, T, 4>) + £(0, <F). (7.10) 

The dual problem is stated as: given Q £ (5 _1 (fl)) 2 , find £ £ V such that 


4($,C) = (Q,$) V$£ V. 


(7.11) 


It can be observed that if is an isolated solution of ( |7.3[ ), then (7.9) and (7.11) are well posed and 
satisfy the a priori bounds 


|e| 2 <c||G|, I©! 


<r c 

II2 — w III III ’ Hl w lll2+a — ^ III'“'III cxllva III Sill 2 — ^ III S* III-1 J HIS Ill2+a — ^ lll^lll-l 5 

where a is the index of elliptic regularity. The discrete linearized problem is defined as: find 0^ G Vh 
such that 

A(0 h , *) = (<?,$) V$£V ;i (7.13) 

where 

A h (G h , $) := 4 h (0 h , $) + B h (% G h , $) + B h (Q h , 4-, $) + £ h (0 h , $). (7.14) 

With this background, Theorem |3.9[ Lemma pi~T| and Theorems |4.2||4~7l can be modihed for the new 
formulation, leading to the applicability of the analysis to a more general form of the von Karman 
equations. We will sketch the proofs of the important results. 


and 


h<c 


ll 2+Q < c 


(7.12) 
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Theorem 7.1. (Well posedness of discrete linearized problem) If'It is an isolated solution of (7.3), 
then for sufficiently small h, the discrete linearized problem (7.13) is well-posed. 


Outline of the proof. Following the proof of Theorem 3.9 we easily arrive at (3.26) using (7.8). To 
estimate |||£7ft.©ft,||| 1 in this case, choose Q = —A EhOh and $ = EhOh in ( 7.11| ) and use (7.13) to obtain 

III-E/i©/Ji = ME h @h - e h , C) + M®h, C - n„c) + B h (% E h Q h , c - n^o + B h (E h o h , it, c - n fc <) 
+ B h (It, E h e h - o h , n^c) + B h (E h Q h -o h , It, n h c) + (G, U h C) 

+ {<L h {E h © - o h , 0 + € h (0 h , C - n h 0) • 


The last term can be estimated using (7.8), Lemmas 3.2 and 3.3 


as 


Kh(E h o - o h , 0 + £ h (Oh, c - n h c)| < Ch |||C||| 2 . 

The remaining terms are estimated as in Theorem |3.9| and result follows. 


(7.15) 

□ 


The next lemma follows as in Lemma 4.1 using (7.8) and hence the proof is skipped. 


Lemma 7.2. (Nonsingularity of perturbed bilinear form) LetTlhlt be the interpolation of It as defined 
in Lemma \3.S\ Then, for sufficiently small h, the perturbed bilinear form defined by 


Mo, $) = A h (0, $) + B h (n h it, o, $) + B h {0, , $) + <t(e, $) 


(7.16) 


is nonsingular on Vi, x Vi, if (7.14) is nonsingular on Vh xVh- 


Theorem 7.3. (Mapping of ball to ball) For a sufficiently small choice of h, there exists a positive 
constant R(h) such that for any O € Vh, 

III© - B h It\\ 2th < R(h) =* MO) - n h $| 2ih < R(h). 


That is, /i maps the ballM^h)(LlhIt) to itself. 


Outline of t he proof. Proceeding as in the proof of Theorem 

/? 


Lemma 7.2 there exists $ £ Vh such that |||$||| 9 h = 1 and 


4.2 


using nonsingularity of Ah(-, ■) and 


HIM©) - n h *| 2 h < AhMO) - n h v, $) 


4 in r~ v “ / -"'-1112 ,h 

= L h {I> - E h It ) + (A h {It, E h i>) - A h (n.hlt, $)) + {B h (It, It, E h I>) - B h (U h H, IL h V, $)) 

5 

+ B h (n h it -o,o- n h it, i>) + (c h (^, E h it) - e h (n h * , $)) =: ^ r<. 


The terms Ti to X4 can be estimated as in the proof of Theorem |4.2[ The last term X5 is estimated 
using (|7.8[), Lemmas 3.2 and|3.3|as: 


|(t h ($,£7 h $)-£ h (n / ,*,$)| < E h $ - $)| + |£ h (tf - IL h It, i)| < Ch |'F||| 2 . (7.17) 

The remaining proof follows exactly same as the proof of Theorem |4.2[ □ 


The existence of solution 'F/ l of (7.5) follows using Theorem 7.3 and satisfies the estimate 


\\It h -IIhIt% th <Ch a . (7.18) 

A contraction result similar to Theorem |4.4| also holds true in this case. The energy estimate follows 
exactly as in the proof of Theorem |4.5| 
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Theorem 7.4. (H 1 estimate) Let 4' and 4/^ be the solutions of (7.3) and (7.5) respectively. Assume 
that 4/ is an isolated solution. Then, for sufficiently small h, it holds 


I® - *fcll 1>h < Ch 2a , 


(7.19) 


where a £ (|, 1] is the index of elliptic regularity. 
Outline of the proof. A use of triangle inequality yields 


II* - */.!U < III* - <h + |||n^ 4 > - * h h h < || 4 / - n ft 4r| + ip _ EhP \L h + \\Ehp\h , (7.20) 


where p = Ilk'S — 4V A choice of Q = — A Ehp and 4> = Ehp in the dual problem (7.11) leads to 


i'VE h p,VE h p) = Ah{E h p , C) = Ah{E h p - pX) + Ah{p, C) 

= Ah(EhP — p, C) + , Ehp — p, C) + Bh{EhP — p , 4/, (f) + ChfEhP — p, C) 

+ AhiUh'S - v, o + A h (* - 4 < h , c - iec) + A h (y,n h c - 0 + MC - n h c) 

+ (B h (% U h Hf - 4v 0 + Bh (- 4 f h , 4', C) - B h { 4 -, 4', 0 + B h (* h) iEO) 
+ (^(iift* - 4-, 0 + <^(4v n fc c - 0) ( 7 . 21 ) 


Combining all the terms related to €h and using (7.8), (7.18) and Lemmas 3.2 3.3 we obtain the 
estimate 

<L h {E h p -pX) + <£h(n„* - 4 -, c) + n„c - C) < Ch 1+a |||C||| 2+Q . (7.22) 


Estimating the remaining terms as in the proof of Theorem |4.6[ the result follows. 


□ 


The Newton’s iterates in this case are defined by 


Akin, $)+B h ( 4-r 1 ,4 >£, 4/”- 1 , $)+£„(*;*, *) = ^(4-r 1 ,4-r 1 , *)+£/,(*) 

The quadratic convergence result follows by a similar proof as in Theorem |4.7[ 


V$eV ft . 

(7.23) 


7.1 Example 3 


In this example, we perform numerical experiments for the problem (7.1)-(7.2) with p/D = 10, over a 
unit square domain. Choose the right hand side load functions such that the exact solution is given 

by 


u{x, y ) = x 2 (l — x) 2 y 2 (l — y) 2 , v(x , y) = sin 2 (7ra;) sin 2 (7ry). 

We consider the same initial triangulation and its uniform refinement process as in Example |5.1[ 
Tables [5] and [6] show the errors and experimental convergence rates for the variables Uh and Vh- 
The computational order of convergences in broken H 2 , H 1 norms are quasi-optimal and verify the 
theoretical results. Also, the order of convergence with respect to L 2 norm is sub-optimal justifying 
the results in m ■ 
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# unknowns 

| U - U h 1 2 , ft . 

Order 

I u - u h \ 1>h 

Order 

||w ^ U h \\ L 2 

Order 

25 

0.101724E-0 

- 

0.129574E-1 

- 

0.469669E-2 

- 

113 

0.391714E-1 

1.3767 

0.275863E-2 

2.2317 

0.957470E-3 

2.2943 

481 

0.195023E-1 

1.0061 

0.767382E-3 

1.8459 

0.252196E-3 

1.9246 

1985 

0.974844E-2 

1.0004 

0.198544E-3 

1.9504 

0.641987E-4 

1.9739 

8065 

0.487399E-2 

1.0000 

0.500990E-4 

1.9866 

0.161298E-4 

1.9928 

32513 

0.243697E-2 

1.0000 

0.125546E-4 

1.9965 

0.403763E-5 

1.9981 


Table 5: Errors and Convergence rates of Uh in broken if 2 , if 1 and L 2 norms 


# unknowns 

\v - V h \ 2 ,h 

Order 

\v - v h h th 

Order 

\\V - V h \\ L 2 

Order 

25 

19.245650 

- 

2.140609E-0 

- 

0.770875E-0 

- 

113 

9.5043692 

1.0178 

0.569978E-0 

1.9090 

0.177898E-0 

2.1154 

481 

5.0549208 

0.9109 

0.161737E-0 

1.8172 

0.482777E-1 

1.8816 

1985 

2.5758938 

0.9726 

0.421546E-1 

1.9398 

0.123930E-1 

1.9618 

8065 

1.2944929 

0.9926 

0.106618E-1 

1.9832 

0.312076E-2 

1.9895 

32513 

0.6480848 

0.9981 

0.267351E-2 

1.9956 

0.781643E-3 

1.9973 


Table 6: Errors and Convergence rates of Vh in broken if 2 , if 1 and L 2 norms 





